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STOIMARY 



The object of this investigation is to study the 
eflect of variation of moisture content in air upon 
the formation of a condensation shock. An ejector 
was used to draw air through a tv;o dimensional con- 
verging-diverging nozzle with glass walls. Moisture 
content of the air was varied a dehumidifier. 
Schlieren pictures of the flow were taken from v/hich 
measurements of the position of the shock were made. 
Also pressure measurements Virere taken along the nozzle. 
The humidity of the inlet air was measured v/ith wet 
and dry bulb thermometers. 

The results obtained by correlation of the data 
for a nozzle of constant angle of divergence are : 

(1) The pressure rise due to a condensation shock 
is a direct function of the specific humidity of the 
inlet air. (Figure VIII). 

(2) For a given inlet temperature, the distance 
of the condensation shock downstream of the throat 
varies inversely as the specific humidity of the air. 
(Figure IV). Prom this relation, it is postulated that 
the position of the area ratio at v/hlch the shock would 
occur can be plotted against specific humidity for a 
given inlet temperature. In this v/ay, the curve will 



serve for nozzles of various contours. (Fic-ure IX). 

(3) The Q'-iount of water vapor condensing out in 
the condensation shock esn be predicted if the inlet 
specific hunidity is known (Fiyure X) . 

(4) The difference between the temperature at 
Y/hich the water vapor becomes saturated under equili- 
brulm conditions and the temperature at which the 
condensation shock occurs, as the humid air nasses 
through the nozzle, is essentially?- a constant of about 
110° p. 

(5) The drop size of the water condensing out 

in the shock is essentially a constant. This is shown 
on Figure XIIo 



INTRODUCTION 



7/hen humid air ie e:^panded by passao;e throuph a 
simple converg"en t-di verpent nozzle, ’the water vaoor 
remains a vapor at the saturation temperature and below 
and continues to expand into the sunersaturated region. 

Part of the water vapor in this meta- stable state con- 

¥ 

denses out rapidly at a certain point along the nozzle. 

At the point of condensation a pressure rise occurs in 
the stream. This phenomena is known as the condensation 
shock. 

A more complete understanding of the condensa- 
tion shock is necessary at the present time due to 
recent increase of flow speeds into the supersonic 
region. Difficulty has been experienced in the design 
of sunersonic wind tunnels and in gas turbine research 
due to the oresonce of condensation shocks. 

Very little information on the condensation shock 
is available in the literature at present. J. I. Yellott 
and C. K. Holland (l)y- (2f' have Investigated the conden- 
sation shock in steam. R. Hermann (3) in 1934-36 studied 
the condensation shock in air. These experiments were 
primarily designed for the provision of data in connec- 
tion with the construction of a workable supersonic wind 
tunnel. Hermann verified previous experiments by showing 

(1) (2) Re'^er to numbered re*'orences in the I^lblior'rL ''dr 
rage 64 . 



that there ia a unique position of the shock in the 
nozzle depending upon the relative humidity of the 
air. K. Oswatitsch (4) published in 1P42 a step by 
step method of computing the pressure as a function 
of distance along a nozzle in the region of the con- 
densation shock. This method was applied to measured 
data of Yellott and Holland and A.k. Binnie and k.W. 
Woods (5), and his own experimental data, and excellent 
correlation was obtained. As mentioned, however, the 

j 

experimental data in this field is small, and the pur- 
pose of our study is to verify previous results and 
obtain additional information on the laws governing the 



condensation shock. 



PROCEDUR.^ AI'ID DESCRIPTION OP AI PAPATUS 



The Schlieren apparatus is an optical -”et>-ioC of 
indicating density gradients of the flow in a noz'^le. 
Th: s a’- oaratus Js -idcquately described in reference 
(6), pp. 26-28 and pp . 38-38. In order to stoo am' 
Dosciblo -noveTient of the shock, the Edgerton }Iigh 
Speed Flash Unit, reference (6), pn. 59 -^- 2 , v/ns used. 
The indicated flash duration of one -five millionths 
of a second .is not attained due to stra:^ wire canaci- 
tance, but the unit provides sufficient speed for the 
purposes of this study. 

A steam e. lector was used to maintain a vacuum of 
about 27 inches of Hg. 'on the nozzle exhaust. This 
insured supersonic flow in the nozzle. A schematic 

vT 

diagram of the apparatus is shown in Figure I. Atmos- 
pheric air may be drawn through valve B and dehumidi- 
fied air through valve A. By adjustment of valves A 

and B, the relative humidity of the air can be varied 

0 

about 405^. The maximum pressure drop through the de- 
humidlfier was 1.5 cm. of Hg. 



Humidity measurement was made by means of a wet 
and dry bulb thermometer at the nozzle inlet. Paired, 
rapid reading thermometers were used. 

As shown in Figure II, reference lines were placed 



at th,G throat and at a station / ,2 inches be lev/ tbo 
throat. This enabled us to accurately ^'leasure the 
position of the shoch and also to deterrnine the fore- 
shortening oi the photograih taken. iressure taps 
v;ere placed along the nozzle at points indicated. 
Pressures v/ere measured by merc.ury manometers. An 
additional pressure tap was located in the 2 ” nozzle 
inlet pipe. 

Simultaneous readings of pressures and of ;vet and 
dry bulb thermometers were made after a steady state of 
flow v/as obtained for the desired relative humidity. 
Either a picture was taken or the nosition of the 
shock measured directly from the camera screen. 



o 




FIGURE I 



DIAGRAMMATIC SKETCH OF APPARATUS 
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table of C^^v'HOLS 



A 

*^Pa 

^Pwv 

S 

H 

ha 

hi 

h\w 

h"' 

J 

k 



m 



a 



- Area, ft. 

- Specific heat at constant pressure of air, 

0o24 R.T.U./lb. Op 

_ Specific heat at constant pressure of v;ater 
vapor, at low pressure, 0.14 B.T.U./lb. °E 

- Acceleration of gravity, 32.2 ft. /sec. ^ 

- Enthalpy of a mixture of air and water vapor 
(plus ice or v;ater if appropriate) B.T.U./lb. 

- Enthalpy of air, B.T.U./lb. 

- Enthalpy of ice, B.T.U./lb. 

- Enthalpy of water vapor, B.T.U./lb. 

- Heat of vaporization or sublimation, P.T.U./lb. 

- Mechanical equivalent of heat, 778 ft. Ibs/P.T.U. 

- Isentropic exponent for the mixture, taken as 1.4 

- Mach No. 

- Lb s . of air 



- Lbs. of water vapor 



m 

P 

P. 






P 



3 00 



sr 



Lbs. of water or ice condensing out in the 

condensation shock 

Pressure of the mixture, lbs. /in. ^ 

Partial pressure of the water vapor, lbs. /in, ^ 

Partial pressure of the air, lbs. /in. ^ 

Saturation pressure for water vapor in thermal 
equilibrium with drops of infinite radius, Ihs./l 

Saturation pressure for water vaoor in thermal 
equilibrium with drops cf radius r, lbs. /in. ^ 

Radius of water droplets, ft. 



t 



Temperature, op 
T - Temperature, absolute 



Tgat “ Teniperature at v/hich the water vapor the 
air becomes saturated, considerinr/; a stab3e 
state exists, as huGild air flows through the 
nozzle, absolute 

V “ Vejocit;/, ft ./sec. 



V 



V 



i 



w or CO 



X 



X’' 



Specific volume of the mixture, ft .3 /lb. 

« 

Specific volume of water at the temperature 
of a drop, ft .5 /lb . 

Specific humidity, lbs. of v/ator* vapor per 
lb . of air 

Actual distance from throat zo condensation 
shock 

Distance from throat to condensation shock 
as measured on a picture 



cr - Surface tension of water, lbs. /ft. 



4 > 



Relative humidity 



Subscripts 



No subscript - Any position in the nozzle 
0 - State at nozzle entrance 



1 

2 



a 



State at an infinitesimal distance 
before the shock 

State at an infinitesimal distance 

after the shock 

Air 



wv 



Water vapor 
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RESULTS 



The data of tables (see Appendix) was 

obtained in sequence. Table III gives the results of 
two sets of runs taken on two different days and is 
arranged' in order of decreasing relative humidity. 

I 

High humidity runs 1, 2, 3, and 4 v/ere obtained using 
steam to humidify the air. The pictures r these 
give qualitative results but the data is considered 
questionable. The pictures obtained during these runs 
follow and are arranged in order of decreasing specific 
humidity. 

Upon reviewing the above data, it was considered 
desirable to obtain more pressure data in the vicinity 
of the shock. Additional pressure taps were added (see 
Figure II) . Table IV presents the results obtained. 
Pressure mxeasurements v/ere taken but no additional rdc- 
tures. Pressures at stations 8 and 4 v/ere unreliable 
due to oblique shocks originating at their pressure 
taps . 

Upon making the plot of x/t (defined on figures 

III and IV) Vs sneclfic humddity, shown in F’igure IV, 

it was believed that x/ 1 v/as a straight line function 

of specific humidity for a given inlet temperature. 

\ 

The data of Table V presents additional inform.ation 
for a lower range of inlet temperatures than had been 
obtained previousl;/. fictures were not taken but 
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I I C TURKS 



Trie ric cures of r*unr, 1, P, 3 rici S are taKen v/1 In 
u'lrk fielT and the condensa t on shoe.; '‘pceai-a l'i'~hl. 
all otlier cic cures are taken ;/i th j.iaht fie 11 and the 
condensation shock shows up as a dark bc-nd. vS'^'cam 
was used to incr'^ase the hunidj ty of runs 1, , " 

and 4. x" is trie distance from the throat to^ the 
ei-i nsa t ion shock ts '"easured cn the picture. ■*<' 

distance fron- the throat to th''' condensation shock ■= 

X = x^ ( 1 .w7) . 






17 . 





Run']., 4 > - 100^, x'"' = 0.34 in. 





3 ^ 



} 



w = 0.0171, x'"‘ - 



Run 4- , - 7 



0.39 in. 






Q O Q 



\ A 




Run 3, ^ -w = 0.0170, x’" = 0.35 in. 



■M 4 * 






Run S, 6 = 56^0, w = 0.01439, x'"' = 0.49 In. 






ir. 





ffun 





0.01265, x'"' 



0,55 in. 




= 0.01190, X'"' = 



Run lo, (j) - 45. G<, 



w 



0 o 60 in. 







0.60 in 



Run 16, 



(/> = ^ 2 %, 



vV 



0.01092, x" = 




Run 10, (^ - 51^,' w = 0.01070, x'" - 0.57 in. 




Run 20, <j) ^ 59 *. w = 0.01027, x'"' = 



0.65 in 






17 




Run 17, </> = 42^:, v/ = 0.01011, x'"' = 0.65 in. 




Run 18, (p = 47<, w = 0.C1C04, x'" = 0.66 in. 




0.59 in. 








Run 23,</>= 34^, w = 0.00964, x'" 



0.75 in. 




Run 21, (j) = olty vv r 0.00951, x'' 



0.c5 in. 







Run 24, ^ = 30. 5<, w = 0.00913, x'"‘ = 0o76 in. 




Run 25, 0 = 27.5'^, w = 0.00856, x'"' = 0.80 in. 




Run 26, 0= 22^, w = 0.00741, x'"‘ = 0.88 Ino 





sho."/D 1 -'lot i. ' ilm n^'’'o:ilcsT ' j 

x/ 1 Vs rolutivo iii.’.r,iil i t,y . This shcv/s ihe s.ir.c 

trend as a sl’-il' r lot of R. f^er'aann (.'■■)• Rela*:n'c 
humidity is not sn absolute quantity. Therefcr-^ ih's 
Diet rennesents qualitative r^'sultg onlv. 

It v/ns decided to use snecJfic l'mmllDt 7 ’' as a 
r arame te:! r s-lacf this is abscli;tc quantity. Sr'.at- 
fic humidity j w, is defined as 11 :)g. of water Vcoor pe 
lb. of ail’. Figure IV shows a plot of x/t Vs a* 
humidity. The inlet temreraturo o^ the air for ;;acr. 
Doint is indicated on the fiauro. It aaoears from 
Fifyjre iv that x/t is a straight line function 
for a given irlet tc-.mperature . 
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DatP for ’'lollin';' curv^r. ] to F on ^Ju-ui-oa V, VT, and VfJ 



was obtained ns ’ndic.'^od ^’n t h--- -^ollowinr tsble. 

TABLE I 



Curve Run Table QQ c O^ y.. (Pp-Pi)/Pn X S y mbol 



1 


1 


II 


.00581 


.00581 


— ±. 1/,/ ^ O 

0.044 


'i 

to 

CO 

< 

o 


0 


2 


2 


II 


.00714 


.00714 


0.051 


0.81 


A 


3 


4 


II 


.00997 








□ 




19 


I 


.00975 


.00985 


0.062 


0.70 


O 




22 


I 


.00982 








A 


4 


5 


II 


.01082 


.01087 


0.073 


0.65 


□ 


» 


16 


I 


.01092 








O 


5 


6 


II 


.01254 


.01259 


0.092 


0.59 


X 




12 


I 


.01265 








^7 


6 


8 


II 


.01358 


.01343 


0.094 


■ 0.55 


V 




9 


I 


.01328 








t 


7 


9 


II 


.01384 


.013814 


Co082 


0.515 


y 


8 


5 


I 


.01588 


.01588 


0.106 


0.43 


o 



As shown In the table, runs with about the same 
specific humidity and Inlet temoerature were grouped, 
where possible, to correlate the data. 



These curves v/ere plotted to represent an instanta- 
neous pressure rise at the position of the shock as de- 
fined below: 

(1) The pressure ratios P/Fq > from the throat down 
to the location of the shock were plotted against dis- 
tance along the nozzle. 

(2) Pressure ratios downstream of the shock v/ere 
plotted. 

(3) These curves were extrapolated to the location 
of the shock as determined from the photographs taken. 

(4) The vertical distance between the two curves 
v/as taken to represent an equivalent pressure rise v/hich 
could be used to represent the condensation shock in 
calculations. This is the value (P2-Pl)/Po* 

Curve VIII shows P2~Pl Vs specific humidity. This 
curve shows that the pressure rise in a nozzle is a di- 
rect function of specific humidity. 

Before proceeding with the discussion of the re- 
sults, the limitation on accuracy of the data presented 
must be mentioned. Steady state conditions for any one 
run were difficult to obtain, but relative humidity never 
varied over Oo5'^. This variation resulted frcm the rise 
of inlet temperature caused b^^ the dehydrating action 
of the silica-gel. Because the pressure measurement 
showed no variation over the duration of a run, it is 



t 



Zl. 




• felt that the specific humidity r mained essentially 
constant . 

kT 

Quick reading, large scale, paired thermometers 
were used in humidity measurements. Air passed over 
wet and dry bulbs at 3,000 ft ./min. which is well 
over the 600 ft. /min. minimum velocity required. 
Relative humidity was obtained by means of reference 
(9) . 

Specific humidity v;as read from Table VII or 
Chart Xtl of reference (8) . The accuracy of the 
recorded specific humidity was about t 0.00CC31bs. 

’water vapor/lb. air. 

Measurements of x, the distance of the shock 
downstream of the throat, were measured to the center 
of the shock to an accuracy of about t 0.01 inch. 

These distances had to be corrected for foreshortening. 

I 

The two reference lines placed 5«2 inches apart on the 

nozzle measured 2.99 inches on the pictures, giving a 

correction factor of 5 . 2 or 1.07 which has to be 

2.99 

applied to measurements taken from pictures. 

Due to the physical nature of the condensi^tion 
shock which occurs over a relatively broad area, the 
measurements were made to the point of maxim.um density 
as mentioned above. It may be shown by comparisons of 
pictures A and B how sensitivity in adjustment of the 



c 

( 



Si'ReAf^ Dir^stTr 



Schlieren Apparatus can vary the apparent width of the 
condensation shock. The condensation shock may be re- 
presented as shown belov/. 




Low Sensitivity Picture A, 
Shock of apparent v/idth I 



Hi^h Sensitivity Picture B, 
Shock of apparent width II 



X 




Picture k - 0.0168, l; 0 \v sensitivity adjustment 




Picture B, u) = 0,01888, Oi e;h sens'f’vity adbistm.ent 





The pressure mea 3 urc."ients obtajned are considered 
accurate to ± 0.5 mm. Hg. Oblique shocks, observed on 
the pictures, from Stations 8 and 4 were believed to 
have resulted In unreliable readings and the pressure 
readings at these stations were not used. 

Upon measuring the areas of the nozzle contour, 
it was found that the measured areas did not agree with 
the nozzle design. Nozzle area measurements were made 
by micrometer to an accuracy of .001 inch. The results 
of these measurements are listed on Table VI in the 
Appendix. A plot of area and area ratio Vs distance 
along the nozzle (Figure XIII) follows Table VI. 



DISCUSSION OP RESULTS 



As brought out in the RI-SULTS, the followinp; 
correlations were obtained for the test nozzle. 

(1) For a given relative humidity, the location 
of the shock can be predicted v/ithin narrow limits. 

(2) For a given inlet temperature and specific 
humidity, the location of the shock can be predicted. 

(5) For a given specific humidity, the pressure 
rise across the shock can be predicted. 

In order that the information of Figure IV may 
be of use for nozzles of various designs. Figure IX 
has been prepared. This figure shows area ratio Vs 
specific humidity, where area ratio is the ratio of 
the area at the position of the shock to the nozzle 
throat area. It is proposed that the results of 
Figure IX will hold for nozzles of various contours. 

The results and conclusions stated to this point 
have been derived directly from the experimental data. 
It is now proposed to investigate the mechanism of the 
condensation shock by analytical methods. In particu- 
lar, the following conclusions will be brought out: 

(1) For a given specific humidity, the amount 
of moisture condensing out at the position of the shoe 



will be •'»redicted. 
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( 2 ) For any specific hinridlty and inlet tempera- 
ture, it win be shown that the difference between the 
temperature at which the air becomes saturated under 
equilibrium conditions and the temperature at which 
condensation occurs is essentially c constant. (Tsat-T 

( 3 ) The order of magnitude of the drop si^e 
formed after the condensation shock 7/111 be ’predicted. 






TIF-OH’'.TTCAL ANALYSIS 



IJathema t ical analvsis of the flov; stream across 
the condensation shock can be made by the use of four 
equations . 

1. Energ:f 

H] + Vi^ - Ho + Vo2 (I) 

2gJ 



2. Continuity 



V 1 - V q 
^ 



since ' 



^2 



( 2 ) 






Home n turn 

r-1 + V.2 



^2 ^2 



2 



V 



IS 



Vog 



or 

Vo = V]_ + gv^ (P. 

irr' 






A 



Svv2 = (Tp, r) 



(3) 



(^) 

( 5 ) 



The use o^" the fourth equation makes two analyti- 
cal methods of approach possible. 

it ff 

(1) Either the drop size r can b-' assim.-d infinite, 
hence fixing the pressure rise due to the condensation 
shock, or 

(2) The observed value of the pressure rls'?. across 



the shock may be used s-nd tht v; Iv-e ol‘ th''- dj-ow 
evaluated. 



. 1 '-e 



Since the second method of approach most nearly 

represents actual conditions, this .';ill be follov/ed. 

>• 

In addition, in ords^r that ■ ctual co .ditioris at the 
shock may be represented as closely as possible, the 
observed value of f^/Po position of ttie shock 

will be used in preference to the theoretical value. 



The energy equation (I) may bo v;ritten: 

mo (han - ha^) + ra „^h ^ - m h „ - m'^h . ^ 

a , 1 2 wvi wvi WV2 \\rv2 i2 



mo + 



■%7Vl 





Sitnplif^ning the left hand side of the equation using: 

^ “ *^vvl ' -%/V2 

^pa ^"^l ~ '^2^ ^'vnri ~ %nj2 ~ ^^^Hvvl ~ ^^~wv2^^'^i2 

m,, + nL,,vi 

V ^ - V ^ 

^2 ^1 

2gJ 



Re-writinc, 



®pa ^/vi ~ ^^2 ^ ~ ^ r72 ~ ^l2^ 

nia ■*■ n^v/vi 




V 



o 



2rJ 



^12 ” ^.vv2 
h = i'i\vv2 



h'^‘ 

hi2 



Substituting above 

^ pa^'^1 ” *^ 2 ^ ~ ^\ rv2 ^ ^ 

*% ‘^rvi 



“ '^’Vv’2 ^ 




V 



2 

1 



2gJ 



“a =na 



■^2^ ^’Vvi ^^,vvi “ ^vv2^ + !^ (^ ) 

ma + 



- ''’z - 

2gJ 



»>a <=na <'^1 " Tg) + (Tj - Tg) t h" 



^a ^vv 1 



O O 

_ V - _ Y- '(7) 
_ . 2 ^ 1 ^ ^ 



So Iv I nr, for fn"', and sort in,' s 1 lb . 



m 



.(1 + V/) (V/'-Y/i + 0.il4(T^-r^) + ;v (0.41) (To-Tj 



iis£. 



Sa''nnle calcul tion ol m'‘ from known nrosruro rise :■ 
shock. 

' Given: 

^o - Psi 

w : lb. /lb. 

. To - Oh 

rp - 1^ - (measured) 



JC 



P 

P 

P 

I 



so 



in. (measured) 
(measured) 



r 1 X P 

1 o 



m 

-o 



(?i)^ 

To ^ 



k-1 



w t 



o 



. 622 (l -f- 



V/ 



) 



.622 



nsi 

OR 



DSi 



To find F 31 . 

Tsi - T 3 Q X - psi 

T , 



(B) 



icroos 



From perfect 



rclatloii usinc; k = 1.4 and irriov/inp- 







Vi = 




f t/ 3 e c 



= ' cr.ft/lb 

?! X 144 

Using eqpation (4) ■' 

V 2 z Vi - gVQ^ (?2 - P]_) X 144 z ft/sec 

^2 ~ ^2 '^1 - 

^1 

The effect of the specific volume of the water ^’anor is 
not included, as the correction is small. 

'^2 - ^2 ^2 ^ 144 - 

R 

h" obtained from Table 5, Reference (10). 

Prom equation (8) 
ra" = Ib/lb 



It is of interest to note that regardless of whether 
h" reorr'sents heat of .fusion or of sublimation, the re ^ 
is litt.le effect on the computed value of m" . Vve chose'' 



*■ c it 



uo : la .e the calcuiation usin~ >'oth values o*^ h'''", it 
will be rointed out loter that it is difficult to ’re- 
dict v/hother ice or \/ater is formed in th^^ shock. 

fhis calculation was carried out ^r the eirnt 
curves on Fio;ures V, VI, and VII. The details cf this 
calculation are shown in the Ap-endix in Table VII. 

The first values of m ' , v/g, and Fgg were obtained usinn 
h equal to the heat of sublimation. For the second 

set of values, h equal to the heat of vaporization 
was used. 

For this computation based on the knov/n pressure 
rise across the shock it was found that the ratio 
^2 Increased as decreased. W 2 is defined as wn - r-'^. 
This information a n shovin on Fl.qure X. It is seen from 
this 4.iqure that w'ith a known inlet specific hur^iditv, 
the lbs. of w'ater vapor per lb air aft^^r the shock can 
be predicted. The m"' used in the calculation of wg for 
Figure X was based upon h’" beJnp, the heat suhlimatlon 




fit.- 

• ^ ■ 



» 



V • 



I 




'I ■> 



It has been found that for steam, v/hen exnandiiip; 
isentropically across the saturation line, that con- 
densation does not occur until the Vi/ilson ne is 
reached. 

This line is located approximately 60 Rtu below 
the saturation line on the Mollier diagram. A simi- 
lar relation v;as sought for humid air expanding isen- 
tropically. 



For a given v/ and Tq, 

^30 - ^^^o 

0.622 



Selecting a pressure ratio F/Pq, 





From the perfect gas relation the temperature of the 
mixture corresr-onding to ^ is 

Fo 



^ ^ U ' L' ' 

^ o 

Prom the calculated T and reference (10), the. saturation 

*? 

pressure, Pgat^ equilibrium conditions .can be com- 

puted. A cut and try process in which various values 
of F/Fq assumed is then used to make F^r ?sat liv- 

ing the temperature at v,hich the air becomes saturated. 
The results of such a computation are shov/n on Figure 
XI. Using competed values of 1^, the following r suits 



v/ere obtained 



3 



TABLE II 



"^sat 



.00581 


496 


.00714 


501 


.00982 


510 


.01082 


514 


.01254 


520 


.01358 


522 


.01384 


523 


.01588 


528 



h 


hat - 


384 


112 


388 


113 


400 


110 


405 


109 


411 


109 


412 


110 


416 


107 


425 


103 



The results Indicate that the condensation shock 
occurs at HO°F. There anoears to be a 

slight trend for the temperature difference to decrease 
v/ith increasing specific humidity. 



4 




lafjQCT 



ISM; 






i'* +-f-{ 






aW3:£ 






-* 

t:;:^ 

^ r 

■ cv 


i-Ltt, 

Lit:: 


4 






"J?"' 




45 . 



We define 


wg - v/2_ 


W 

- m 


(9) 


then, v/g 


- 0.622 


P32 


(10) 






P2 




or. P32 = 


WgPgX ■ 


* psi 


(11) 


0.622 







Values of Pgg were computed for the eight runs 
for which m'“* v/as computed. It v;as found upon comparing 
Pgg computed from formula (11), that Pgg much greater 

than the saturation pressure for equilibrium conditions 
corresponding to the temperature Tg. 

Now v/e will review the events occuring in the flow 
of the humid air through the nozzle. At the nozzle in- 
let the water vapor in the air is in a supa^ heated 
state. The mixture of water vapor and air is expanded 
isentropically through the nozzle and at a temperature 
of about 110 degrees below the temperature at v/hich 
condensation should have occured, some of the v/ater 
vapor condenses out. As mentioned above, the vapor 
pressure after the shock is greater than the equili- 
brium vapor pressure at that temperature. The vapor 
i'^ressure under stable conditions is that corresponding 
to vapor in equilibrium v/lth drops of infinite radius. 

The Von-Helmholtz relation shows that if the pressure 



•^•Reference (ll) , pages 434 - 437 



of the vv .ter vanor Is rreater than the sa^-uration 
pressure under equilibrui.Ti conditions, th-:n the v/ater 
vapor Is in equillbruim v/ith drons of less than in- 
finite radius. 

^ By means of the Von-He Imholtz equation, (see 
Aopendlx pp. 62), Figure XII, which follows, was drawn, 
ligure XII shows, for various drop sizes, the pressure 
of the v/ater vapor in equilibruim with those drops for 
a given temperature. Calculated points from the eight 
runs v/hich were calculated in detail (Table VII) are 
shown on Figure XII. These points are plotted using 
T 2 and Pgg as arguments. The circled points correspond 
to Pg 2 computed on the basis tha t orms and the "x" 
points correspond to Pg^ computed on the basis that^^Bl 
forms. The four low specific .humidity points lie at a 
pressure below the triple point pressure of 0.0888 psi. 
Under stable conditions, ice would be formed in the shock’ 
at these lov/ pressures. The vapor pressures of the four 
high humidity points lie above the triple point pressure 
and probably water droplets would form in the shock. 
However, it is difficult to predict in this metastable 
statq whether ice or water droplets form. For the four 
high humidity runs, it appears that a drop size of about 
3 X 10 ^ feet is formed in the shock. 



It is possible that 
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FIGURE XB 

PLOT SHOVilNO THE SIZE OF WATER DROPLETS 




PRESSURE FOR WATER VAPOR IN EOUl LI B R U IN VVITW 
P.R0P5, OF RADIUS ^HOy^ NI - PSI , 



4C. 



u I'le t. L ’ tub le stat of rfubcooled ./a ter exists bolov; the 
trinle point pressure. If this is true, it arrears from 
Figure XII that for all snoci bic hu-nidi ties the drop size 
for’-;iln ■' in the shock is of a radius of about 3 x 10“^ feeto 
Because of our inability to predict the exact nature of 
this phenomena, both the circled and ’’x" points are shown. 

It must be mentioned that in the comoutation of 
drop size it was found that the various reliable authori- 
ties ^ive an appreciable variation in thv values of 
Pg^ at teraperatures below 32®F. These discrepencies 
increase with decreasing temperatures. Also values of 
sui'face tension are not definitely known at these low 
temperatures. 

In conclusion, it is obvious that further exploration 
of this sub.iect is necessary in order to present a cornnlete 
nicture. However, it is honed that the results of our 
experiments w.1 ] 1 contribute to the knov/ledee of this 
sub iec t . 



RIi^CO.ilv'ENDATIONS 



(1) In preference to the conventional tv/o dimen- 
sional nozzle, it would be better to use a half nozzle 
contour with a sliding straight wall section fitted 
with a pressure tap and a micrometer screw for accurat 
longitudinal measurement of the position at v-hich the 
pressure is measured. 

(2) liore accurate humidity measurements are de- 
sirable. Chemical methods can give any desired degree 
of accuracy. 

(.■ 5 ) It is desirable to cover a wider range of 
inlet temperatures and specific humidities. 

(4) Verify the nroposed correlation that a .plot 
of area ratio against specific humidity for a given 
inlet temperature will determine the position of the 
condensation shock for nozzles of various contours. 
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VI 
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61 
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62 



librium v/ith drops of various 
radii . 



TABLE III 



Run 


1 


2 


3 


4 


5 


UPITS 


% Rel.Hurao 


100 


36.8 


77 


75.5 


65 




D.B. 


90.8 


82.9 


80.4 


80.8 


82.7 


Op 


V/.B. 


90.8 


79o6 


74.9 


74.7 


73.4 


Op 


Po 


75.7 


75o70 


75.7 


75.65 


75.65 


cm. Hg . 


^10 


42o55 


40o70 


40.45 


40.25 


40.20 


tt 


Pg 


36o95 


35.90 


35.35 


35.45 


33.85 


n 


^8 


32.55 


32 . 30 


32 . 30 


32.40 


32.25 


TI 


P 7 


25o20 


24.30 


23.90 


23.75 


23.55 


n 


^2 


12.80 


12.00 


11.70 


11.60 


11.35 


n 


CaJ 










.01588 


lb 3. W. V./ 

lbs. air 


■ Picture 


Yes 


Yes 


Yes 


Yes 


"'G. 




Pio/Po 










.532 




V^o 










.448 




V^o 










.427 




P 7 /P 0 










,312 




V^o 










.150 














o41 




-p 











.84 





TABLEHr (Conte «d) 



Run 


6 


7 


8 


9 


10 


11 


% R.H. 


I 63 


61 


55 


51 


51 


50 


D.B. 


78 


83.2 


84 0 3 


84.9 


78o2 


78.3 


W.B. 


69 


73.0 


72.4 


71.4 


65.8 


65.5 


Po 


76.15 


75o65 


75.6 


75.6 


76.15 


76.15 


PlO 


40.20 


40.15 


40.15 


40.05 


40.20 


40 . 10 


P9 


35.20 


33.45 


32.60 


52 . 30 


52.30 


32o25 


Ps 


30.95 


32.25 


32 . 35 


32.25 


30.85 


30 . 90 


P7 


22.95 


23.50 


23.20 


22.95 


22.10 


22.05 


Ps 


11.15 


11.30 


11.20 


11.10 


10.70 


10.65 


CJ 


.01313 


.01513 


.01439 


.01328 


.01070 


.01053 


Picture 


Yes 


Yes 


Yes 


Yes 


Yes 




Pic/Po 


.528 


.532 


.532 


.530 


.528 


.526 


P9/P0 


.436 


o442 


.431 


.427 


.424 


o424 


Ps/Po 


.407 


.427 


.428 


.427 


.405 


o406 


P7/P0 


.302 


.311 


.307 


.304 


.290 


.290 


P2/P0 


.146 


.1496 


.148 


.147 


.141 


.141 


x" 


o47 


.43 


.49 


.53 


.57 




47 1 


o96 


.88 


1.00 


1.08 


I0I6 





TABLSin (Cont. »d) 



Run 


12 


13 ' 


14 


15 


16 


17 


% R.H. 


47 


46 


44 


43.6 


42 


42 


D.B. 


85.7 


79 


80 0 5 


36 


84.6 


82.2 


'U . B . 


70.8 


64.8 


65.7 


59 . 8 


68 . 2 


66.2 


Po 


75.40 


76.05 


75.95 


75.30 


75.25 


75.75 


PlO 


39.95 


40.05 


40.00 


39.85 


39.95 


39.95 


^9 


32.15 


32.20 


32.15 


32.00 


32.00 


32.05 


*’8 


32.30 


30.30 


30o25 


31 oOO 


30.40 


28.60 


Py 


22.60 


21.85 


21.80 


22.05 

i 


22.00 


21.55 


P 2 


10.90 


10.55 


10o55 


10 06^5 


10 . 60 


10.45 


00 


.01265 


.00992 


.00998 


0 OII 9 O 


.01092 


.01011 


Picture 


Yes 


Yes 


No 


Yes 


Yes 


, Yes 


Pio/Po 


o530 


.527 


o527 


.530 


.531 


.527 


P9/P0 


o426 


.424 


o423 


o425 


.425 


.423 


Pe/Po 


.428 


.398 


.398 


c411 


.404 


.378 


P7/P0 


.300 


.287 


.287 


c292 


.292 


.284 


Pg/Po 


.145 


.139 


.139 


.1415 


.141 


.138 


-k'- 

X 


.55 


.59 




. 60 


.60 


.65 


cr 


1.12 


1.20 




lo22 


1.22 


1.33 




TABLEIU (Pont, <d) 



10 19 20 21 



nun 
io R.H. 


xo 

40 


39.5 


39 


37 


35 


34 


D.B. 


83.5 


83.0 


34.9 


34 0 2 


86 


87.2 


W.B. 


6 0 • 5 


66.2 


67.1 


66 .0 


66.8 


66.9 


Po 


75.75 


75.15 


75.75 


75.75 


75.65 


75.65 


^10 


39.90 


39.75 


39.90 


39.90 


39.90 


59.75 


Pg 


32.05 


31.95 


32.05 


32.05 


32.05 


31.90 


^8 


28.40 


30.05 


25. . 60 


28.60 


27 . 90 


25 o 45 




21.60 


21.65 


21.65 


21.65 


21.75 


21.70 


^2 


10.45 


10.40 


10.50 


10 0 50 


10 0 50 


10.40 


oO 


.01004 


,00975 


0OIO27 


.00951 


.00982 


.00954 


Picture 


Yes 


No 


Yes 


Yes 


Yes 


Yes 


Pio/Po 


.527 


. 530 


.527 


.527 


.526 


.526 


P9/P0 


.423 


o 425 


o 423 


.423 


.424 


.422 


Ps/Po 


.378 


o 400 


.378 


.378 


.369 




P7/P0 


.285 


.288 


.286 


.286 


.288 


.287 


P2/P0 


.138 


.1385 


.139 


.139 


.139 


.133 


x“ 


.66 




.65 


.65 


o 65 


c 75 


^ t 


1.35 




lo 33 


1.33 


1.33 


1.53 



TABLEIg (Cont. 



24 25 26 



% R.H. 


30.5 


27.5 


/CO 

22 


D.B. 


89.7 


90.2 


92.4 


W.B. 


66 . 5 


66.7 


66.5 


^0 


75.25 


75.20 


75.0 


PlO 


39 .65 


39.65 


39o50 


Pg 


31.80 


31.80 


31.65 


Ps 


25.30 


24.95 


24.45 


P? 


21.85 


21.85 


21.40 


Ps 


10 . 35 


10.35 


10.25 


Ol> 


.00913 


.00856 


.00741 


Picture 


Yea 


Yes 


Yes 


^lo/Po 


.527 


.527 


o526 


P9/P0 


.423 


.423 


.422' 


Pa/Po 


.336 


.332 


. 326 


P7/P0 


.290 


.290 


.286 


Ps/Po 


.138 


.138 


.136 


x" 


.76 


oSO 


.Sn 


x/ 1 


1.55 


1.63 


1.8 



TABLE IV 



Hun 1 2 3 4 5 UlIITS 



% Re 1. Ham. 


24.0 


27.5 


34.5 


37 . 5 


41.0 


Op 


D.B. 


82.4 


04.5 


85.4 


85 . 5 


85.0 


Op 


VV.B. 


59.8 


62.9 


65.8 


67.0 


68. C 


Op 


Po 


74o95 


75.05 


75.40 


7 5.60 


75.70 


cm. Hg. 


PlO 


40.00 


40.1 


40.25 


40.40 


40 . 50 


tt 




37.95 


38.00 


3Bo20 


38.30 


38.40 


It 


Pg 


32.70 


32.70 


32.85 


32 . 95 


33.00 


II 




27o20 


27 . 30‘ 


27.50 


27.90 


28.65 


M 


Ps^ 


31.00 


51.25 


32.40 


33.60 


34.60 


II 


P 7— 


22o85 


23.95 


25.15 


2600 


26.50 


tt 


P 7 


21.10 


21.80 


21.75 


22.00 


22.25 


ft 


Pe 


16.80 


17 . 30 


18.15 


18.70 


18.06 


It 


P 4 * 


14.30 


14.65 


15.10 


15.40 


15.06 


It 


Pg 


9.90 


10.15 


10 . 35 


10.55 


10.70 


tt 


Pio/Po 


.554 


.535 


. 534 


.535 


• 5o5 




Pgj/Po 


.505 


.506 


.v506 


.507 

1' 


.507 




P9/P0 


.436 


.436 


.435 


.436 


.436 




Psi/Po 


.363 


.364 


.365 


.369 


.378 




Pyj/Po 


o305 


.318 


.333 


.344 


.350 




P7/P0 


.282 


.290 


.289 


.291 


.294 




Ps/Po 


.224 


.231 


o'241 


.248 


.246 




P2/P0 


.120 


.135 


ol37 


.140 


.141 




OL) 


c00581 


,00714 


.00923 


.00997 


.01082 


lbs. '.V. V./ 
lbs. air 



V a lue s no t 



ured 



TABLE IV (Cont.’d) 



Run 


6 


7 


8 


9 


10 


% R.H. 


45.0 


49.0 


55.5 


58.0 


98.0 


D.B. 


33.5 


83 . 6 


82.8 


32.0 


87.0 


W.B. 


68.5 


69.4 


70.6 


70. B 


86. 5 


Po 


75.80 


75.90 


76.10 


76.10 


76.10 


^10 


40.60 


40.65 


40.85 


40 .70 


41.85 




38.50 


38.55 


38.55 


38.60 


40.55 


P 9 


33.10 


33.15 


33.35 


33.60 


37 . 30 




30.20 


31.10 


32.90 


34.60 


34.60 


Pg* 


35.80 


36.55 


39.30 


40.50 


38.30 


P7I 


26.90 


27.30 


27o05 


27.09 


28.80 


P? 


22.65 


22.80 


23.25 


23.30 


25.10 


Ps 


18.70 


18.80 


19.05 


19.15 


21.30 


P 4 '> 


15.85 


15.95 


I 60 O 5 


16.20 


17.60 


P2 


10.90 


10.95 


11.10 


11.10 


12.70 


Pio/Po 


.535 


.536 


.537 


.535 


c550 


P94/P0 


.508 


c508 


o507 


.507 


.534 


P9/P0 


.436 


.436 


.430 


.441 


.490 


P8*/Po 


.398 


o410 


o432 


.455 


o455 


Pv^Po 


.354 


.360 


.356 


.355 


.379 


P7/P0 


.299 


.300 


.306 


.306 


.330 


Pe/Po 


.246 


.248 


.250 


.252 


o2B0 


P2/P0 


.144 


.14^ 


cl46 


.146 


ol67 


CjJ 


0OI254 


.01233 


.01358 


o01384 


.02764 






7aliu'. s not need 




TABLE. 



Run 


ic Re 1, Hum. 


<jJ 


D.B. 


V/.B. 


x' 




1 


63o0 


.0107 


72.4 


53.8 


.49 


1.00 


2 


57 


.0099 


73.4 


63.2 


.53 


1.08 


5 


50 


.0088 


73.5 


61.4 


.57 


1.16 


4 


46 


e0083 


74.2 


61.0 


.59 


1.20 


5 


45 


.0030 


74.0 


60.4 


.60 


1.225 


6 


42.5 


.0076 


74.5 


60.0 


.61 


lo245 


7 


41 


.0075 


74.6 


59.8 


.62 


1.265 


8 


36.5 


.0070 


75.7 


59.4 


.65 


1 . 327 


Q 

V/ 


. 51 


c00645 


76.6 


59.0 


.69 


1.410 


10 


30 


.00585 


77.6 


58.5 


.73 


lo490 


11 


30 


.00610 


78.2 


59.0 


.74 


1.510 


12 


28.5 


.00600 


78.7 


59.0 


.79 


lo610 



TnHLL VI 



DaT/. of nozzle AREiiS AND AREA R.iTIOS 



Stationc aro s’oaced inchea rirart ( F1 rn.irn 

follo'win,,^) . The nozzle is 0.50 inches thick. 



Stati on Area , In.^ A//vtp 



2 


.-H2 


1.804 


3 


.■1135 


1 . 590 


4 


.384 


1.567 


5 


. 356 


1.453 


6 


. 328 


1.310 


7 


.2975 


1.214 


8 


.271 


1.107 


9 


.248 


1.012 


10 (Throat ) 


.2455 


1.000 


11 


.2155 


] .002 


12 


.246 


1.004 


13 


.2V9 


lol39 


14 


.335 .. 


1.37 



1 
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• 






- ■ " 














CURVE NO 


1 


Z 


5 


4 


5 


6 


7 


8 




ft 


lARS 


75.03 


7565 


75.7 
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The calcrlation of the vapor pressure of water 
vapor in thermal equilibrium v/ith drops of various 
radii for various temperatures: 



Ln Pgj, = 2 cr v^ 

"rRT 



'Ahere R - 85.8 and (T is 
dynes/cm. rriven by: 



the surface 

kT 



tension in 



or = 75.64 - 0.1391t - cOOOSt'^ (t in °C) 

^^300 ^3 given by the equation belov;; 

logio 21.075 - 290v5.59 - 4.71754 TogiQT 

T 

'Ahere T is expressed in ‘^Kelvin. 

^ SCO psi. 
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